coexistence diameter in two-dimensional colloid-polymer mixtures 
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We demonstrate that the law of the rectilinear coexistence diameter in two-dimensional (2D) 
mixtures of non-spherical colloids and non-adsorbing polymers is violated. Upon approach of the 
critical point, the diameter shows logarithmic singular behavior governed by a term tint, with t 
the relative distance from the critical point. No sign of a term t^^ could be detected, with f3 the 
critical exponent of the order parameter, indicating a very weak or absent Yang- Yang anomaly. Our 
analysis thus reveals that non-spherical particle shape alone is not sufficient for the formation of a 
pronounced Yang- Yang anomaly in the critical behavior of fluids. 
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Critical phenomena have been studied extensively for 
decades. Nevertheless, a few important issues remain 
controversial. One longstanding problem concerns the 
critical behavior of the so-called coexistence diameter 
of fluids undergoing a phase transition from liquid to 
gas (the coexistence diameter is defined as the average 
density of the two coexisting phases). In 1964, Yang 
and Yang showed that the divergence of the constant- 
volume specific heat at the critical point implies that 
either (fip/dT"^ or d^^/dT^ or both diverge Here, 
T is the temperature, p the pressure, and ji the chem- 
ical potential. A remarkable consequence, realized only 
recently by Fisher and co-workers is that the diver- 
gence of (Pji/dT^ implies that the coexistence diameter 
5 gains an additional term t^^ which, assuming Ising uni- 
versality, dominates the previously recognized term 
I^Q. As usual, t is the relative distance from the critical 
point; a and (3 are the critical exponents of the specific 
heat and the order parameter, respectively. The coexis- 
tence diameter near the critical point thus reads as 



5a#0 = Per (1 + A2pt'^^ + Ai^^t^-" + A^t) 



(1) 



with Per the critical concentration, and non-universal am- 
plitudes Ai. The divergence of (fp/dT'^ is called a Yang- 
Yang (YY) anomaly, and 7^ if one is present. 
For symmetric fluids, such as the Widom-Rowlinson mix- 
ture d^p/dT^ remains finite, in which case A213 ~ 
0. However, realistic fluids are typically asymmetric, in 
which case a YY- anomaly cannot be ruled out. 

Indeed, evidence for a YY-anomaly was found in simu- 
lations of asymmetric three-dimensional (3D) fluids 0,13 ■ 
Experiments on propane and carbon dioxide also show 
evidence for a weak YY-anomaly Q, while for He-3 the 
situation is less clear However, for 3D fluids, the 

analysis is extremely difficult. Assumiiig 3D Ising univer- 
sality with a w 0.109 and /3 w 0.326 [llj, there are two 
singular terms in Eq.(^ which need to be distinguished 
not only from each other, but also from the leading an- 
alytic background term Ait. A direct observation of the 
term t"^^ is consequently very difflcult. Only very re- 



cently, after it was recognized that the amplitudes Ai^a 
and Ai are coupled, could both singular terms in Eq.QJ 
be resolved from experimental data [l^ . 

Nevertheless, due to the competition between terms 
in Eq.l^, investigations of the coexistence diameter in 
3D fluids remain challenging. An attractive alternative, 
where the competition between singular terms is less se- 
vere, is to consider a fluid where the specific heat diverges 
logarithmically at the critical point, implying a = 0. The 
critical behavior of the diameter is then given by 



[1 + A2i3t^^ + Aq tint + Alt) 



(2) 



In other words, the term is replaced by a (much 

weaker) logarithmic singularity |l3l |. In order to detect 
the YY-anomaly, one thus needs to distinguish power law 
behavior from logarithmic behavior, which is expected to 
yield a more pronounced numerical signature. To this 
end, one could consider a 2D fluid, where a — and 
(3 = 1/8 (assuming 2D Ising universality). Note that 
mean-field systems do not qualify, despite having a = 0. 
For mean-field systems, a = corresponds to a finite 
discontinuity in the specific heat, yielding a purely rec- 
tilinear diameter |l3| . Still, even though the numerical 
analysis of a 2D fluid with Ising critical behavior may 
become more simple, this is no guarantee that a YY- 
anomaly will be found. The absence or presence of a 
YY-anomaly is a non-universal feature: A2/3 may well 
be zero! It is not completely clear which features in a 
fluid determine the strength of the YY-anomaly. Obvi- 
ously, the fluid must be asymmetric. In addition, there 
are experimental indications that molecular shape and 
symmetry, in particular departures from spherical form, 
are important contributing factors 0. From these con- 
siderations, it appears that the "minimal" fluid in which 
a YY-anomaly may most easily be found should (1) be 
two-dimensional and exhibit a 2D Ising critical point (2) 
be asymmetric, and (3) contain non-spherical particles. 

In this work, we will investigate a fluid with pre- 
cisely these properties, and focus on the critical behav- 
ior of its diameter. We consider a 2D version of the 
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colloid-polymer model of Asakura and Oosawa (AO) [l4| . 
but generalized to non-spherical colloids. In the orig- 
inal AO model, colloids and polymers are treated as 
spheres in 3D, assuming hard-core interactions between 
colloid-colloid and colloid-polymer pairs, while polymer- 
polymer pairs can interpenetrate freely (the AO model is 
thus clearly asymmetric). Since the polymers may over- 
lap freely, their translational entropy is increased signifi- 
cantly when the colloids group together. Hence, there is 
an effective (depletion) attraction between the colloids. 
Provided polymer concentration and size are sufficiently 
large, the attraction is strong enough to drive phase sep- 
aration in the AO model, whereby the system splits-up 
into a colloid-rich (polymer-poor) phase, and a colloid- 
poor (polymer-rich) phase. As expected for systems with 
short-ranged interactions, the corresponding unmixing 
critical point belongs to the 3D Ising universality class 
[l5| . Hence, it is anticipated, although one should plan 
to check, that the AO model in 2D will exhibit a 2D Ising 
critical point. The effect of non-spherical particle shape 
is incorporated by modeling the colloids not as spheres, 
but as line segments. Our minimal model is thus a 2D 
mixture of colloidal line segments of length 1^, and effec- 
tive polymer disks of diameter Cp, interacting via AO po- 
tentials. In other words, overlaps between line segments, 
as well as overlaps between line segments and polymer 
disks are forbidden, while the polymer disks may overlap 
freely. 

The aim of this work is to check if a YY-anomaly in this 
model can be found. We will do so using computer simu- 
lation and finite size scaling (FSS) in the grand canonical 
ensemble. In this ensemble, the total area of the system 
A, the temperature T, and the colloid (polymer) fugac- 
ity Zc (zp) are fixed, while the number of colloids {Nc) 
and polymers (iVp) fluctuates. The thermal wavelength 
is set to unity, such that Zp reflects the average concen- 
tration Np/A a pure phase of polymers would have (recall 
that such a phase is simply an ideal gas). The remain- 
ing lengths are expressed in units of ^c- The colloid-to- 
polymer size ratio is set to ap/lc = 0.95. At the coex- 
istence colloid fugacity, it is expected that the mixture 
phase separates into a colloid poor phase (the gas) and 
a colloid rich phase (the liquid), provided the polymer 
fugacity exceeds a critical value Zp^cr- The phase sepa- 
ration is thus driven by Zp, which therefore plays a role 
analogous to inverse temperature in gas-liquid transitions 
of simple fluids. The relative distance from the critical 
point is written as i = Zp/zp,ci — 1, and pQ (pl) denotes 
the concentration pc = Nc/A of colloids in the gas (liq- 
uid) phase. A natural order parameter is the density gap 
A = (pl ~ Pg)/2, while the coexistence diameter reads 
as S = (pl + Pg)/2. In the limit t — > 0, we expect critical 
power law behavior A oc for the order parameter, and 
Eq.(|2Jl for the coexistence diameter. 

For our minimal model, the critical point must be lo- 
cated first, in particular the value of Zp.cr- Next, 2D Ising 




FIG. 1: upper frame: Binodals of the 2D colloid-polymer mix- 
ture of this work, obtained using mean-field theory (MF) and 
simulation (SIM); circles mark critical points, lower frame. ■ 
Order parameter A as function of t, obtained using simulation 



and FSS; a value Zn 



■ 3.877 in t was used. 



universality of the critical point must be established, af- 
ter which the diameter can be investigated. To "guide" 
the simulations, the phase diagram is obtained approxi- 
mately first, using a simple mean-field (MF) theory based 
on a free- volume approach for 2D rod-polymer mixtures. 
The free- volume fraction and the free energy for a pure 
system of 2D lines are derived from scaled particle the- 
ory The free energy is exactly the same as the one 
obtained within Onsager's second virial theory jl^. The 
resulting binodal is shown in Fig. ^ The theory predicts 
the critical point at p, 
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nematic ordering of the colloids is predicted in the direct 
vicinity of the critical point: the coexisting phases are 
isotropic. Since the theory ignores critical fiuctuations, it 
is expected that underestimates the true value Zp^cr 
significantly. Nevertheless, the theoretical result is im- 
portant because it provides an indication in which regime 
the (time consuming) simulations need to be carried out. 
The simulations are performed in the grand canonical 
ensemble, on a 2D square of size A = L x L using peri- 
odic boundary conditions. We measure the distribution 
Pl(Pc|zp, Zc), defined as the probability of observing a 
system with colloid concentration pc , at fugacities Zp and 
Zc, with L the system size. The insertion and removal of 
particles is performed using a cluster move |15j, com- 
bined with a biased sampling scheme [l7| to overcome 
the free energy barrier separating the phases, and his- 
togram reweighting ll^. To obtain a single distribution, 
around 100 CPU hours for a small system (L = 22), and 
350 h for a large system {L = 30), are required. 

A standard route to obtain Zp^cr in simulations, is 
to measure the L-dependence of the cumulant U4 = 
{vn?) / {rn^) along some path in the (pc, Zp)-plane. Here, 
m — Pc ^ {pc), and (■) denote grand canonical averages. 
The cumulant becomes system-size independent at the 
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FIG. 2: Cumulant analysis near the critical point. The top 
frame shows U4 as function of Zp for various system sizes L 
as indicated. The lower frame shows the cumulant slope li 
at Zp^ci as function of L. 



0.00 0.02 0.04 0.06 ^ 0.08 0.10 0.12 



order parameter scaling curve 


_ y=1+x/2 


0.15 






coexistence diameter scaling curve 




y' 






0.10 






0.05 




y'-x72 








0.00 



0.00 O.02 0.04 X' 0.06 0.08 0.10 

FIG. 3: upper frame: Order parameter scaling curve y = f{x) 
(solid line). Following convention the scaling curve is 
raised to a negative exponent with <j!> — l/P, where the 2D 
Ising value /3 = 1/8 was used. Also shown is the small x 
limiting form y = 1 + x/2 (dashed line), lower frame: Co- 
existence diameter scaling curve y' — g{x') (solid line). The 
dashed line shows the small x' limiting form y' = x' /2. 



critical point [l9j. Plots of U4 as function of Zp for dif- 
ferent system sizes L are expected to show a common 
intersection point, leading to an estimate of Zp,cr- More- 
over, the value of the cumulant at the intersection 
point is universal, and this gives an indication of the uni- 
versality class. The result is shown in the upper frame of 
Fig. |2] The data were obtained using the colloid fugacity 
at which (to^) is maximized 1^. From the intersections, 
we obtain Zp ^r — 3.881 ± 0.005, where the error reflects 
the scatter in the various intersection points. At the in- 
tersection point w 0.85, which is very close to the 
accepted 2D Ising value ~ 0.856 fi^l (horizontal Hue 
in Fig. 121). While this already suggests 2D Ising universal- 
ity, additional confirmation is obtained from the critical 
exponents. We consider (3 and ly, with j/ the critical expo- 
nent of the correlation length. Here, v is extracted from 
the cumulant slope Yi = dUi/dzp at the critical value 



of Zp. One expects that Yi oc L^^'^, with L the system 
size. The lower frame of Fig.[21shows Yi as function of L, 
where the above estimate of Zp.cr was used. The line is 
a linear fit through the origin, which describes the data 
very well, and thus confirms the 2D Ising value — I. 
To obtain /3, we apply the FSS algorithm of Ref.Q, using 
system sizes L = 20 — 30. Starting with Zp significantly 
above its critical value, the algorithm proceeds by plot- 
ting U4 as function of the average colloid concentration 
(pc)- The resulting plot reveals two minima, located at 
p~ and p"*", with respective values Q~ and at the 
minima. Defining the quantities Qmin = {Q~^ + Q^)/2, 
X = <3minln(4/e(3min), and y = (p+ - p-)/(2A), the 
points (a;, y) obtained for different system sizes L should, 
in the limit far away from the critical point, collapse onto 
the line y — 1 + x/2. Recall that A is the order parame- 
ter in the thermodynamic limit at the considered fugacity 
Zp, precisely the quantity of interest, which may thus be 
obtained by fitting until the best collapse onto 1 -I- a::/2 
occurs. In the next step, Zp is chosen closer to the criti- 
cal point, the points (x, y) are calculated as before, but 
this time A is chosen such that the new data set joins 
smoothly with the previous one, yielding an estimate of 
the order parameter at the new fugacity. This procedure 
is repeated as closely as possible to the critical point, 
where A vanishes. The output of the algorithm, A as 
function of Zp, is then fitted to A oc i'', in order to esti- 
mate Zp^cr and (3. We obtain Zp^cr = 3.877 ±0.001, which 
is consistent with the (less precise) analysis of Fig. |5] 
Shown in the lower frame of Fig. ^ is the order parame- 
ter as function oft, on double logarithmic scales. The line 
has slope (3 = 1/8, and confirms the 2D Ising exponent in 
the simulation data. The FSS algorithm also yields the 
order parameter scaling curve y = f(x) 0, with x and 
y defined as above, shown for completeness in the upper 
frame of Fig. |21 Away from the critical point [x — ^ 0), 
the scaling curve has the limiting form y 1 + x/2\ at 
the critical point, the scaling curve diverges. The sig- 
nificance of the scaling curve is its universal character: 
all systems with a 2D Ising critical point should yield a 
scaling curve for the order parameter similar to the one 
shown here 0- 

At this point, sufficient evidence for 2D Ising univer- 
sality has been provided. Since the order parameter is a 
scalar, and since the interactions are short-ranged, there 
are in any case no theoretical grounds to contemplate a 
different universality class. Hence, we will now consider 
the coexistence diameter. To obtain the diameter, the 
FSS algorithm of Ref. is used. Similarly, also for the 
coexistence diameter, a scaling curve y' — g{x') is gen- 
erated, based on different quantities x' and y' defined in 
Ref.H In contrast to the order parameter, the scaling 
curve of the diameter is not universal. For our model, the 
corresponding curve is shown in the lower frame of Fig.O 
For small x' , it correctly approaches the exact limiting 
form y' = a;'/2 [tI. The curvature at x' ^ already sug- 
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FIG. 4: upper frame: Coexistence diameter S as function of 
t (circles). The dashed line is a fit to Eq.Q. lower frame: n 
as function of \nt (circles). The straight line (dashed) con- 
firms the logarithmic nature of the divergence. In both plots, 
Zu.cr = 3.877 in t was used. 



gests singular behavior. According to Eq.Q, this may 
reflect the YY-anomaly, or logarithmic behavior, or both. 
To quantify this, the diameter itself is shown in the up- 
per frame of Fig. 01 The quality of the data is such that 
5 can be resolved down to t « 0.0015. A fit to Eq.Q 
yields Per = 0.9270 ± 0.0006, « 0, ^0 = -0.29 ± 0.01 
and Ai = 0.76 ±0.03, where the error reflects the scatter 
resulting from the range over which the fit is performed. 
The fit describes the data perfectly well, without the need 
for a term t^^ ^ indicating logarithmic singular behavior. 
In other words, despite the non-spherical particle shape 
in our model, the present analysis does not reveal a YY- 
anomaly. For completeness, in the lower frame of Fig. ^ 
the derivative n = d6/dt is plotted as function of Int. 
In case of singular behavior, k is expected to diverge as 
t — > 0. The (logarithmic) divergence is clearly visible. 
Finally, by combining the coexistence diameter and or- 
der parameter data, the binodal in the thermodynamic 
limit was constructed, see the top frame of Fig. ^ 

In conclusion, we have shown that 2D colloid-polymer 
mixtures, with non-spherical colloids, do not display a 
pronounced YY-anomaly, at least not for the colloid-to- 
polymer size ratio q = 0.95 considered by us. Although 
the diameter becomes singular upon approach of the crit- 
ical point, the singularity is logarithmic, and well de- 
scribed by the theoretically expected term tint. While 
it has been suggested that non-spherical particle shape 
may be an important contributing factor to the forma- 
tion of the YY-anomaly, this seems not to be the case 
for our 2D model. If a YY-anomaly is present in our 
2D model nevertheless, it is very weak, and negligible 
down to t « 0.0015 accessible in our simulations. In 
contrast, our results may be compatible with the very 
recent Ref. J^, where it is argued that a YY-anomaly is 
expected when p* = pcrO-i is small, where a; represents 
the typical interaction volume. For our model, ai ~ q^, 
and so p* « 0.8. This value even exceeds p"* w 0.75 of the 



Widom-Rowlinson mixture |3] , for which no YY-anomaly 
was detected either . In order to detect a YY-anomaly 
in our 2D model, it seems that smaller size ratios q are re- 
quired; this could be a topic for further simulations. It is 
tempting to speculate if the 2D model considered in this 
work can also be realized experimentally. Colloidal parti- 
cles, due to their mesoscopic size, pose many advantages 
over atomic fluids. This has already enabled the investi- 
gation of critical phenomena in 3D, whereby the particles 
are visualized directly using confocal microscopy. Other 
applications may be found in order-disorder phase transi- 
tions in adsorbed monolayers of atoms or small molecules 
at surfaces. 
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